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Abstract: We study the production of a top-quark pair in association with one and two 
vector bosons, ttV and tiVV with V = 7 , Z, at the LHC. We provide predictions at 
next-to-leading order in QCD for total cross sections and top-quark charge asymmetries as 
well as for differential distributions. A thorough discussion of the residual theoretical uncer¬ 
tainties related to missing higher orders and to parton distribution functions is presented. 
As an application, we calculate the total cross sections for this class of processes (together 
with tiH and titi production) at hadron colliders for energies up to 100 TeV. In addition, 
by matching the NLO calculation to a parton shower, we determine the contribution of tiV 
and tiVV to hnal state signatures (two-photon and two-same-sign-, three- and four-lepton) 
relevant for tiH analyses at the Run II of the LHC. 
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1 Introduction 

With the second run of the LHC at 13 TeV of centre of mass energy, the Standard Model 
(SM) is being probed at the highest energy scale ever reached in collider experiments. At 
these energies, heavy particles and high-multiplicity final states are abundantly produced, 
offering the opportunity to scrutinise the dynamics and the strength of the interactions 
among the heaviest particles discovered so far: the W and Z bosons, the top quark and 
the recently observed scalar boson [1, 2], The possibility of measuring the couplings of the 
top quark with the W and Z bosons and the triple (quadruple) gauge-boson couplings will 
further test the consistency of the SM and in case quantify possible deviations. In addition, 
the couplings of the Higgs with the W and Z bosons and the top quark, which are also 
crucial to fully characterise the scalar sector of the SM, could possibly open a window on 
Beyond-the-Standard-Model (BSM) interactions. 

Besides the study of their interactions, final states involving the heaviest states of 
the SM are an important part of the LHC program, because they naturally lead to high- 
multiplicity final states (with or without missing transverse momentum). This kind of 
signatures are typical in BSM scenarios featuring new heavy states that decay via long 
chains involving, e.g., dark matter candidates. Thus, either as signal or as background 
processes, predictions for this class of SM processes need to be known at the best possible 
accuracy and precision to maximise the sensitivity to deviations from the SM. In other 
words, the size of higher-order corrections and the total theoretical uncertainties have to 
be under control. In the case of future (hadron) colliders, which will typically reach higher 
energies and luminosities, the phenomenological relevance of this kind of processes and the 
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impact of higher-order corrections on the corresponding theoretical predictions are expected 
to become even more relevant [3]. 

In this work we focus on a specific class of high-multiplicity production process in the 
SM, i.e., the associated production of a top-quark pair with either one {ttV) or two gauge 
vector bosons {ttVV). The former includes the processes + tiW~), ttZ and 

tt'y, while the latter counts six different final states, i.e., ttW~^W~, ttZZ, ttyy, 
ttW^Z and ttZ^. In addition, we consider also the associated production of two top-quark 
pairs {tttt), since it will be relevant for the phenomenological analyses that are presented 
in this work. 

The aim of our work is twofold. Firstly, we perform a detailed study of the predictions 
at fixed NLO QCD accuracy for all the ttV and ttVV processes, together with ttH and tttt 
production, within the same calculation framework and using the same input parameters. 
This approach allows to investigate, for the first time, whether either common features 
or substantial differences exist among the theoretical predictions for different final states. 
More specifically, we investigate the impact of NLO QCD corrections on total cross sections 
and differential distributions. We systematically study the residual theoretical uncertainties 
due to missing higher orders by considering the dependence of key observables on different 
definitions of central renormalisation and factorisation scales and on their variations. NLO 
QCD corrections are known for ttH in [4-7], for ttj in [ 8 , 9], for ttZ in [9-13], for ttW^ 
in [9, 13-15] and for ttti in [16, 17]. NLO electroweak and QCD corrections have also 
already been calculated for ttH in [18-20] and for ttW^ and ttZ in [20]. Moreover, in 
the case of ttH, NLO QCD corrections have been matched to parton showers [21, 22] and 
calculated for off-shell top (anti)quarks with leptonic decays in [23]. In the case of ti'j, NLO 
QCD corrections have been matched to parton showers in [24]. For the ttVV processes a 
detailed study of NLO QCD corrections has been performed only for ttjj [25, 26]. So far, 
only representative results at the level of total cross sections have been presented for the 
remaining ttVV processes [3, 17]. When possible, i.e. for ttV, ttH and ttyy, our results have 
been checked against those available in the literature in previous works [9, 13, 14, 17, 20- 
22, 24, 25], and we have found perfect agreement with them. This cross-check can also be 
interpreted as a further verification of the correctness of both the results in the literature and 
of the automation of the calculation of NLO QCD corrections in MadGraph5_aMC@NLO. 

Secondly, we perform a complete analysis, at NLO QCD accuracy including the match¬ 
ing to parton shower and decays, in a realistic experimental setup, for both signal and 
background processes involved in the searches for ttH at the LHC. Specifically, we con¬ 
sider the cases where the Higgs boson decays either into two photons {H —)• 77 ), or into 
leptons (via H —?■ WW*, H —?■ ZZ*, H —?• which have already been analysed by 

the CMS and ATLAS collaborations at the LHC with 7 and 8 TeV [27-29]. In the first 
case, the process ttyy is the main irreducible background. In the second case, the processes 
tiW+W-, tiZZ, tiW^Z are part of the background, although their rates are very small, 
as we will see. However, ttW~^W~ production, e.g, has already been taken into account 
at LO in the analyses of the CMS collaboration at 7 and 8 TeV, see for instance [27]. A 
contribution of similar size can originate also from tttt production [30], which consequently 
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has to be included for a correct estimation of the background.^ Furthermore, depending 
on the exact final state signature, the ttV processes can give the dominant contribution, 
which is typically one order of magnitude larger than in tiVV and tttt production. 

In this work, the calculation of the NLO QCD corrections and the corresponding event 
generation has been performed in the MadGraph5_aMC@NLO framework [17]. This code 
allows the automatic calculation of tree-level amplitudes, subtraction terms and their inte¬ 
gration over phase space [31] as well as of loop-amplitudes [9, 32, 33] once the relevant Feyn¬ 
man rules and UV/ R 2 counterterms for a given theory are provided [34-36] . Event genera¬ 
tion is obtained by matching short-distance events to the shower employing the MC@NLO 
method [37], which is implemented for PythiaG [38], PythiaS [39], HERWIG6 [40] and 
HERWIGH— 1-[41]. The reader can find in the text all the inputs and set of instructions 
that are necessary to obtain the results presented here. 

The paper is organised as follows. In section 2 we present a detailed study of the predic¬ 
tions at NLO QCD accuracy for the total cross sections of tiV, ttVV and tttt production. 
We study their dependences on the variation of the factorisation and renormalisation scales. 
Furthermore, we investigate the differences among the use of a fixed scale and two possible 
definitions of dynamical scales. Inclusive and differential iL-factors are also shown. As 
already mentioned above, these processes are backgrounds to the tiH production with the 
Higgs boson decaying into leptons, which is also considered in this work. To this purpose, 
we show also the same kind of results for tiH production. In addition, in the case of tiV 
and tiH, we provide predictions at NLO in QCD for the corresponding top-charge asymme¬ 
tries and in order to investigate the behaviour of the perturbative expansion for some key 
observables, we also compute tiVj and tiHj cross sections at NLO in QCD. Such results 
appear here for the first time. In section 2 we also study the dependence of the total cross 
sections and of global AT-factors for tiVV and tiV processes as well as for tiH and titi 
production on the total energy of the proton-proton system, providing predictions in the 
range from 8 to 100 TeV. 

In section 3 we present results at NLO accuracy for the background and signal relevant 
for tiH production. In subsection 3.1 we consider the signature where the Higgs decays 
into photons. In our analysis we implement a selection and a definition of the signal region 
that are very similar to those of the corresponding CMS study [27]. For the signal and 
background processes ti'j'j, we compare LO, NLO results and LO predictions rescaled by a 
global flat AT-factor for production only, as obtained in section 2. We discuss the range of 
validity and the limitations of the last approximation, which is typically employed in the 
experimental analyses. In subsection 3.2 we present an analysis at NLO in QCD accuracy for 
the searches of tiH production with the Higgs boson subsequently decaying into leptons (via 
vector bosons), on the same lines of subsection 3.1. In this case, we consider different signal 
regions and exclusive hnal states, which can receive contributions from titi production and 
from all the tiV and tiVV processes involving at least a heavy vector boson. Also here, 

^Triple top-quark production, tttW and titj, a process mediated by a weak current, is characterised by 
a cross section that is one order of magnitude smaller than tttt at the LHC and it is usually neglected in 
the analyses. 
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we compare LO, NLO results and LO predictions rescaled by a global flat iT-factor for 
production only. In section 4 we draw our conclusions and present an outlook. 

2 Fixed-order corrections at the production level 

In this section we describe the effects of fixed-order NLO QCD corrections at the production 
level for ttV processes and tiH production (subsection 2.1), for ttVV processes (subsection 
2.2) and then for titi production (subsection 2.3). All the results are shown for 13 TeV 
collisions at the LHC. In subsection 2.4 we provide total cross sections and global iL-factors 
for proton-proton collision energies from 8 to 100 TeV. With the exception of tt’j’j, detailed 
studies at NLO for tiVV processes are presented here for the first time. The other pro¬ 
cesses have already been investigated in previous works, whose references have been listed 
in introduction. Here, we (re-)perform all such calculations within the same framework, 
MadGraph5_aMC@NLO, using a consistent set of input parameters and paying special 
attention to features that are either universally shared or differ among the various processes. 
Moreover, we investigate aspects that have been only partially studied in previous works, 
such as the dependence on (the definition of) the factorisation and renormalisation scales, 
both at integrated and differential level. To this aim we define the variables that will be 
used as renormalisation and factorisation scales. 

Besides a fixed scale, we will in general explore the effect of dynamical scales that 
depend on the transverse masses {mT,i) of the final-state particles. Specifically, we will 
employ the arithmetic mean of the of the final-state particles (/Tq) and the geometric 
mean (/r^), which are defined as 

Ht 1 v—' 

j=i,Ar(+i) 

( V'" 

Aig := mT,i (2.2) 

\i=l,N J 

In these two definitions N is the number of final-state particles at LO and with V(-|-l) in 
eq. (2.1) we understand that, for the real-emission events contributing at NLO, we take 
into account the transverse mass of the emitted parton.^ There are two key aspects in 
the definition of a dynamical scale: the normalisation and the functional form. We have 
chosen a “natural” average normalisation in both cases leading to a value close to rrit when 
the transverse momenta in the Born configuration can be neglected. This is somewhat 
conventional in our approach as the information on what could be considered a good choice 
(barring the limited evidence that a NLO calculation can give for that in first place) can 
be only gathered a posteriori by explicitly evaluating the scale dependence of the results. 
For this reason, in our studies of the total cross section predictions, we vary scales over 

^This cannot be done for /Tg; soft real emission would lead to fig ~ 0. Conversely, fia can also be defined 
excluding the partons from real emission and, in the region where mr.i’s are of the same order, is numerically 
equivalent to fig. We remind that by default in MadGraph5_aMC@NLO the renormalisation and 
factorisation scales are set equal to Ht/2. 
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a quite extended range, /ic /8 < /r < 8 /ic- More elaborate choices of even-by-event scales, 
such as a CKKW-like one [42] where factorisation and renormalisation scales are “local” and 
evaluated by assigning a parton-shower like history to the final state configuration, could be 
also considered. Being ours the first comprehensive study for this class of processes and our 
aim that of gaining a basic understanding of the dynamical features of these processes, we 
focus on the simpler definitions above and leave possible refinements to specific applications. 

All the NLO and LO results have been produced with the MSTW2008 (68% c.l.) PDFs 
[43] respectively at NLO or LO accuracy, in the five-flavour-scheme (5FS) and with the 
associated values of q;<j. tiW~^W~ production, however, has been calculated in the four- 
flavour-scheme (4FS) with 4FS PDFs, since the 5FS introduces intermediate top-quark 
resonances that need to be subtracted and thus unnecessary technical complications. 

The mass of the top quark has been set to mt = 173 GeV and the mass of the Higgs 
boson to niH = 125 GeV, the GKM matrix is considered as diagonal. NLO computations 
are performed by leaving the top quark and the vector bosons stable. In simulations at 
NLO+PS accuracy, they are decayed by employing MadSpin [44, 45] or by PythiaS. If 
not stated otherwise photons are required to have a transverse momentum larger than 
20 GeV {pt{i) > 20 GeV) and Frixione isolation [46] is imposed for jets and additional 
photons, with the technical cut Rq = 0.4. The fine structure constant a is set equal to its 
corresponding value in the G^-scheme for all the processes.^ 

2.1 tiV processes and tiH production 

As first step, we show for tiH production and all the tiV processes the dependence of the 
NLO total cross sections, at 13 TeV, on the variation of the renormalisation and factorisation 
scales pr and pf. This dependence is shown in fig. 1 by keeping p = pr = Pf and varying 
it by a factor eight around the central value p = pg (solid lines), fi, = fi,a (dashed lines) and 
fj, = mt (dotted lines). The scales pa and Hg are respectively defined in eqs. (2.1) and (2.2). 

As typically pa is larger than pg and mt, the bulk of the cross sections originates 
from phase-space regions where as{pa) < ois{pg), as{mt). Gonsequently, such choice gives 
systematically smaller cross sections. On the other hand, the dynamical scale choice fig 
leads to results very close in shape and normalisation to a fixed scale of order mt- 

Driven by the necessity of making a choice, in the following of this section and in the 
analyses of section 3 we will use fig as reference scale. Also, we will independently vary 
fif and fir by a factor of two around the central value fig, fig/2 < fif,fir < “^fig, in order 
to estimate the uncertainty of missing higher orders. This generally includes, e.g., almost 
the same range of values spanned by varying fi = fij. = fif hy a factor of four around the 
central value fi = fia, fia/^ < fi < 4/ia (cf. fig. 1 ) and thus it can be seen as a conservative 
choice. In any case, while certainly justified a priori as well as a posteriori, we stress that 
the fi = fig choice is an operational one, i.e. we do not consider it as our “best guess” 
but just use it as reference for making meaningful comparisons with other possible scale 
definitions and among different processes. 

®This scheme choice for a is particularly suitable for processes involving W bosons [47]. Anyway, in our 
calculation, no renormalisation is involved in the electroweak sector, so results with different values of a 
can be obtained by simply rescaling the numbers listed in this paper. 
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ttH 
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ttw^ 

tty 


NLO 

rt + 6.0% +2.1% 

OZZ.Z_g_4^^ _2.6% 

ofTO ^+10.3% +2.0% 
^''^■'^-11.7% -2.5% 

CAA n + 13.0% +1.7% 
b44.«ii g% 3% 

2746+i4'2% 

+ 1.6% 
-1.9% 

LO 

/17A (:?+35.5% +2.0% 
^/0.0 24 2% -2.1% 

71 n 0“l”36.1% +2.0% 
/1U.0 24 5% -2.1% 

q+ 28.1% +1.7% 
ozo.y 2o.4% -1.8% 

2100^^®'^^“ 
^-LUU 24 5% 

+ 1.8% 
-1.9% 

A-factor 

1.10 

1.23 

1.22 

1.31 



Table 1. NLO and LO cross sections for tiV processes and tiH production at 13 TeV for ^, = fig. 
As already stated in the text, with ttW^ we refer to the sum of the tiW~^ and tiW~ contributions. 
The first uncertainty is given by the scale variation within fig/2 < fif,fir < ‘2fig, the second one by 
PDFs. The relative statistical integration error is equal or smaller than one permille. 

Using the procedure described before, in table 1 we list, for all the processes, LO 
and NLO cross sections together with PDF and scale uncertainties, and iL-factors for the 
central values. The dependence of the LO and NLO cross sections on /r = /r,. = /rj- is also 
shown in fig. 2 in the range fJ-g/8 < fi < 8fig. As expected, for all the processes, the scale 
dependence is strongly reduced from LO to NLO predictions both in the standard interval 
lJig/2 < fi < 2fig as well as in the full range fig/8 < fi < 8//^. For tt'y process (upper- 
plots in figs. 1 and 2), we find that in general the dependence of the cross-section scale 
variation is not strongly affected by the minimum pT of the photon, giving similar results 
for pt{i) > 20GeV and PTil) > SOGeV. As already stated in section 1, with tiW^ we 
refer to the sum of the tiW~^ and ttW~ contributions. 

We now show the impact of NLO QGD corrections on important distributions and we 
discuss their dependence on the scale variation as well as on the definition of the scales. For 
all the processes we analysed the distribution of the invariant mass of the top-quark pair 
and the Pt and the rapidity of the (anti)top quark, of the top-quark pair and of the vector or 
scalar boson. Given the large amount of distributions, we show only representative results. 
All the distributions considered and additional ones can be produced via the public code 
MadGraph5_aMC@NLO. 

For each figure, we display together the same type of distributions for the four different 
processes: tby, ttH, tiW^ and ttZ. Most of the plots for each individual process will be 
displayed in the format described in the following. 

In each plot, the main panel shows the specific distribution at LO (blue) and NLO 
QGD (red) accuracy, with p = pf = pr equal to the reference scale pg. In the first inset 
we display scale and PDF uncertainties normalised to the blue curve, i.e., the LO with 
p = Pg. The mouse-grey band indicates the scale variation at LO in the standard range 
Pg/2 < Pf,Pr < 2,pg, while the dark-grey band shows the PDF uncertainty. The black 
dashed line is the central value of the grey band, thus it is by definition equal to one. The 
solid black line is the NLO QGD differential AT-factor at the scale p = pg, the red band 
around it indicates the scale variation in the standard range Pg/2 < Pf,Pr < 2,pg. The 
additional blue borders show the PDF uncertainty. We stress that in the plots, as well as 
in the tables, scale uncertainties are always obtained by the independent variation of the 
factorisation and renormalisation scales, via the reweighting technique introduced in [48]. 
The second and third insets show the same content of the first inset, but with different 
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Figure 1. Cross sections at 13 TeV. Comparison of the NLO scale dependence in the interval 
/r'^/8 < fj, < 8^'^ for the three different choices of the central value fj,g, "m-t- The upper plot 
refers to tij production, the lower plot to tiW^, ttZ and tiH production. 

scales. In the second panel both LO and NLO have been evalnated with /x = /x^, in the 
third panel with /x = m*. 

The fonrth and the fifth panels show a direct comparison of NLO QCD predictions nsing 
the scale [Xg and, respectively, fia and mt- All cnrves are normalised to the red cnrve in the 
main panel, i.e., the NLO with /x = fig. The monse-grey band now indicates the scale varia¬ 
tion dependence of NLO QCD with fi = fig. Again the dashed black line, the central valne, 
is by dehnition eqnal to one and the dark-grey borders represent the PDF nncertainties. 
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Figure 2. LO and NLO cross sections at 13 TeV. Scale dependence in the interval /r'^/8 < fj, < S/i”^ 
with /i° = fig. The npper plot refers to ti^ prodnction, the lower plot to tiW^, ttZ and tiH 
production. 

The black solid line in the fonrth panel is the ratio of the NLO QCD predictions at the scale 
fia and jjLg. The red band shows the scale dependence of NLO QCD predictions at the scale 
fia) again normalised to the central valne of NLO QCD at the scale fig, denoted as R{fia)- 
Bine bands indicate the PDF nncertainties. The hfth panel, R{mt), is completely analogons 
to the fonrth panel, bnt it compares NLO QCD predictions with fig and mt as central scales. 


We start with hg. 3, which shows the distribntions for the invariant mass of the top- 
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Figure 3. Differential distributions for the invariant mass of top-quark pair, m{ti). The format of 
the plots is described in detail in the text. 

quark pair {m{ti)) for the four production processes. From this distribution it is possible 
to note some features that are in general true for most of the distributions. As can be 
seen in the fourth insets, the use of ^ = /Tq leads to NLO values compatible with, but 
systematically smaller than, those obtained with // = /r^. Conversely, the using ^ = rrit leads 
to scale uncertainties bands that overlap with those obtained with fj, = fig. By comparing 
the first three insets for the different processes, it can be noted that the reduction of the 
scale dependence from LO to NLO results is stronger in ttH production than for the tiV 
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Figure 4. Differential distributions for the Pt of top-quark pair, pT{ti). The format of the plots is 
described in detail in the text. 

processes. As we said, all these features are not peculiar for the m{tt) distribution, and are 
consistent with the total cross section analysis presented before, see hg. 1 and table 1. From 
hg. 3 one can see that the two dynamical scales fig and pa yield flatter iF-factors than those 
from the hxed scale mt, supporting a posteriori such a reference scale. While this feature 
is general, there are important exceptions. This is particular evident for the distributions 
of the Pt of the top-quark pair {pT{ti)) in hg. 4, where the differential iF-factors strongly 
depend on the value of pritt) for both dynamical and hxed scales. The relative size of QCD 
corrections grows with the values of Priti) and this effect is especially large in tiW^ and 
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Figure 5. Representative kinematical configurations for tiW final state. At LO (left) a high-p^ ti 
pair recoils against the W boson. At NLO (right), the dominant configuration is the one where the 
jet takes most of the recoil and the W boson is soft. 

ti'j production. In the following we investigate the origin of these large iC-factors. 

Top-quark pairs with a large Pt originate at LO from the recoil against a hard vector 
or scalar boson. Conversely, at NLO, the largest contribution to this kinetic configuration 
emerges from the recoil of the top-quark pair against a hard jet and a soft scalar or vector 
boson (see the sketches in fig. 5). In particular, the cross section for a top-quark pair with a 
large pr receives large corrections from (anti)quark-gluon initial state, which appears for the 
first time in the NLO QCD corrections. This effect is further enhanced in ttW^ production 
for two different reasons. First, at LO ttW^ production does not originate, unlike the other 
production processes, form the gluon-gluon initial state, which has the largest partonic 
luminosity. Thus, the relative corrections induced by (anti)quark-gluon initial states have 
a larger impact. Second, the emission of a LF collinear to the final-state (anti)quark in 
qg —)• tiW^q' can be approximated as the qg —)• ttq process times a g —?■ q'W^ splitting. For 
the W momentum, the splitting involves a soft and collinear singularity which is regulated 
by the W mass. Thus, once the W momentum is integrated, the qg —>■ tiW^q' process yields 
contributions to the Priti) distributions that are proportional to a* log^ \pT{tt)/mw].^ The 
same effect has been already observed for the Pt distribution of one vector boson in NLO 
QCD and EW corrections to Z and ZZ bosons hadroproduction [49-51]. 

The argument above clarifies the origin of the enhancement at high of the ti pair, yet 
it raises the question of the reliability of the NLO predictions for tiV in this region of the 
phase space. In particular the giant AT-factors and the large scale dependence call for better 
predictions. At first, one could argue that only a complete NNLO calculation for ttV would 
settle this issue. However, since the dominant kinematic configurations (see the sketch on 
the right in fig. 5) feature a hard jet, it is possible to start from the ttVj final state and 
reduce the problem to the computation of NLO corrections to ttVj. Such predictions can 
be automatically obtained within MadGraph5_aMC@NLO. We have therefore computed 

^In ttZ the same argument holds for the q —>■ qZ splitting in qg —>■ tiZq. However, the larger mass of 
the Z boson and especially the presence of the gluon-gluon initial state at LO suppress this effect. 
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Figure 6. Comparison between differential distribution of the tt transverse momentum in tiW^ as 
obtained from calculations performed at different orders in QCD. The blue and red solid histograms 
are obtained from the tiW^ calculation at LO and NLO respectively. The dashed histograms are 
obtained from the tiW^j calculation at LO (light blue, purple, and mouse-grey) and at NLO 
(green), for different minimum cuts (50, 100, 150 GeV) on the jet p-r- The lower inset shows the 
differential i^T-factor as well as the residual uncertainties as given by the tiW^j calculation. 

results for different minimum Pt for the extra jet both at NLO and LO accuracy. In fig. 6 
we summarise the most important features of the tiW^{j) cross section as a function of 
the Pxitt) as obtained from different calculations and orders. Similar results, even though 
less extreme, hold for ttZ and ttH final states and therefore we do not show them for sake 
of brevity. In fig. 6, the solid blue and red curves correspond to the predictions of Pritt) as 
obtained from ttW^ calculation at LO and NLO, respectively. The dashed light blue, purple 
and mouse-grey curves are obtained by calculating tiW^j at LO (yet with NLO PDFs and 
Us and same scale choice in order to consistently compare them with NLO ttW^ results) 
with a minimum cut for the jets of 50, 100, 150 GeV, respectively. The three curves, 
while having a different threshold behaviour, all tend smoothly to the tiW^ prediction at 
NLO at high pT{tt), clearly illustrating the fact that the dominant contributions come from 
kinematic configurations featuring a hard jet, such as those depicted on the right of fig. 5. 
Finally, the dashed green line is the Priti) as obtained from ttW^j at NLO in QCD with 
a minimum p^ cut of the jet of 100 GeV. This prediction for pT{tt) at high p^ is stable and 
reliable, and in particular does not feature any large iC-factor, as can be seen in the lower 
inset which displays the differential iL-factor for tiW^j production with p^ cut of the jet 
of 100 GeV. For large pritt), NLO corrections to tiW^j reduce the scale dependence of LO 
predictions, but do not increase their central value. Consequently, as we do not expect large 
effects from NNLO corrections in ttW^ production at large pT{tt), a simulation of NLO 
ttV+jets merged sample a la FxFx [52] should be sufficient to provide reliable predictions 
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13 TeV a[fb] 

ttHj 


ttZj 

ttWj 


NLO 

1 /I Q q+3.3% 

+3.0% 

-3.6% 

ooA 7+6.6% +2.8% 

^'^'^■'-13.4% -3.2% 

202 

zuz.»_igg% 

+1.4% 

-1.1% 

LO 

1 74 ir+57.s% 
^^^■^-33.9% 

+2.8% 

-2.9% 

04*4 1+68.2% +2.7% 
Z4:0. l _34 o% _2.8% 

197 fi +53 ' 7 % 

+1.5% 

-1.5% 

iL-factor 

0.85 


0.95 

1.03 



Table 2. Cross sections with prij) > 100 GeV. The renormalisation and factorisation scales are 
set to /ig of ttV. The (N)LO cross sections are calculated with (N)LO PDFs, the relative statistical 
integration error is equal or smaller than one permil. 

over the full phase space. 

For completeness, we provide in table 2 the total cross sections at LO and NLO accuracy 
for tiW^j, as well as ttZj and ttHj production, with a cut ptU) > 100 GeV. At variance 
with what has been done in fig. 6, LO cross sections are calculated with LO PDFs and the 
corresponding as done in the rest of the article. 

The mechanism discussed in detail in previous paragraphs is also the source of the giant 
iC-factors for large pritt) in ti'y production, see fig. 4. This process can originate from the 
gluon-gluon initial state at LO, however, the emission of a photon involves soft and collinear 
singularities, which are not regulated by physical masses. When the photon is collinear to 
the final-state (anti)quark, the qg — )• tt'jq process can be approximated as the qg — ?• ttq 
process times a q ^ q'j splitting. Here, soft and collinear divergencies are regulated by 
both the cut on the Pt of the photon (p™*) and the Frixione isolation parameter Rq. We 
checked that, increasing the values of p™* and/or Rq, the size of the iC-factors is reduced. 
It is interesting to note also that corrections in the tail are much larger for p = pg than 
p = Pa- This is due to the fact that the softest photons, which give the largest contributions, 
sizeably reduce the value of pg, whereas pa is by construction larger than 2pT{tt}. This 
also suggests that pg might be an appropriate scale choice for this process only when the 
minimum p^ cut and the isolation on the photon are harder.^ 

In figs. 7 and 8 we show the p^ distributions for the top quark and the vector or scalar 
boson, Prit) and Pt(H), respectively. For these two observables, we find the general features 
which have already been addressed for the m{tt) distributions in fig. 3. 

In fig. 9 we display the distributions for the rapidity of the vector or scalar boson, y{V). 
In the four processes considered here, the vector or scalar boson is radiated in different ways 
at LO. In ttH production, the Higgs boson is never radiated from the initial state. In ttZ 
and tt'j production, in the quark-antiquark channel the vector boson can be emitted from 
the initial and final states, but in the gluon-gluon channel it can be radiated only from 
the final state. In tiW^ production, the W is always emitted from the initial state. The 
initial-state radiation of a vector boson is enhanced in the forward and backward direction, 
i.e., when it is collinear to the beam-pipe axis. Consequently, the vector boson is more 
peripherally distributed in ttW^ production, which involves only initial state radiation, 

®Assuming mT{t) ~ rriTit) and my( 7 ) = , the the ratio pa/pg increases by increasing prit) and, 

when mT{t) > decreases by increasing p^^. Moreover, under the same assumption, pa = Pg at 

rriTit) = Pt’*. For these reasons, especially for large PtRE), pg may underestimate the value of the scale. 
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Figure 7. Differential distributions for the Pt of top-quark, The format of the plots is 

described in detail in the text. 

than in tij and especially tiZ production. In ttH production, large values of |y(I^)| are 
not related to any enhancement and indeed the y{V) distribution is much more central 
than in ttV processes. These features can be quantified by looking, e.g., at the ratio 
•= %(\y\ = ^)/%i\y\ = 3 )- At LO we find, r{W) ~ 5, r{-y) ~ 8.5, r{Z) ~ 17.5 
and r(H) ~ 40. As can be seen in the first three insets of the plots of fig. 9, NLO QCD 
corrections decrease the values of r{V) for ttW^ and tt'y production, i.e. the vector bosons 
are even more peripherally distributed {r{W) 3.5, r(7) 5.5). A similar but milder effect 

is observed also in tiZ production (r(Z) ~ 16). On the contrary, NLO QCD corrections 
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Figure 8. Differential distributions for the pr of the vector or scalar boson, Pt{V). The format of 
the plots is described in detail in the text. 

make the distribution of the rapidity of the Higgs boson even more central {r{H) ~ 53). 
In fig. 9 one can also notice how the reduction of the scale dependence from LO to NLO 
results is much higher in ttH production than in ttV type processes. Furthermore, for this 
observable, iL-factors are in general not flat also with the use of dynamical scales. From a 
phenomenological point of view, this is particularly important for ttW^ and ti^, since the 
cross section originating from the peripheral region is not extremely suppressed, as can be 
seen from the aforementioned values of r(LF) and r('j). 


- 15 - 
































































































y(H) 




Figure 9. Differential distributions for the rapidity of the vector or scalar boson, yiV). The format 
of the plots is described in detail in the text. 

In fig. 10 we show distributions for the rapidities of the top quark and antiquark, y{t) 
and y(t). In this case we use a slightly different format for the plots. In the main panel, 
as in the format of the previous plots, we show LO results in blue and NLO results in red. 
Solid lines correspond to y{t), while dashed lines refer to y(t). In the first and second inset 
we plot the ratio of the y{t) and y{t) distributions respectively at NLO and LO accuracy. 
This ratio is helpful to easily identify which distribution is more central (peripheral) and if 
there is a central asymmetry for the top-quark pair. Also here, although it is not shown in 
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Figure 10. Differential distributions for the rapidity of the top quark and antiquark, y{t) and y{t). 
the plots, itT-factors are not in general flat. 

In the case of ti production the central asymmetry, or the forward-backward asymmetry 
in proton-antiproton collisions, originates from QCD and EW corrections. At NLO, the 
asymmetry arises from the interference of initial- and final-state radiation of neutral vector 
bosons (gluon in QCD corrections, and photons or Z bosons in EW corrections) [53-58]. 
Thus, the real radiation contributions involve, at LO, the processes pp —)• ttZ and pp —)• tty, 
which are studied here both at LO and at NLO accuracy. As can be seen from fig. 10, tty 
production yields an asymmetry already at LO, a feature studied in [59]. The tiZ production 
central asymmetry is also expected to be non vanishing at LO, but the results plotted in 
fig. 10 tell us that the actual value is very small. The asymmetry is instead analytically 
zero in tiW^ (ttH) production, where the interference of initial- and final- state W (Higgs) 
bosons is not possible.® 

At NLO, all the ttV processes and the ttH production have an asymmetry, as can be 

®In principle, when the couplings of light-flavour quarks are considered non-vanishing, initial-state radi¬ 
ation of a Higgs boson is possible and also a very small asymmetry is generated. However, this possibility 
is ignored here. 
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13 TeV Ac [%] 

ttW^ tty 

LO 

-3.93t°J6 ± 0.03 

NLO 

2.90tn7 -om ± 0.07 -1.79t°;™ ± 0.06 

13 TeV Ac [%] 

tiH tiZ 

LO 

-0.12t°:[;[ ± 0.03 

NLO 

1 .00t!]iS] -o;o4 ± 0.02 0.85^°;f7 ± 0.03 


Table 3. NLO and LO central asymmetries for ttV-type processes and tiH production at 13 TeV 
for /r = p,g. The first uncertainty is given by the scale variation within fJ-g/2 < p,f,p,r < 2/ig, the 
second one by PDFs. The assigned error is the absolute statistical integration error. 

seen in fig. 10 from the ratios of the y{t) and y(t) distributions at NLO. In the case of ttW^ 
production the asymmetry, which is generated by NLO QCD corrections, has already been 
studied in detail in [15]. In all the other cases it is analysed for the first time here. 

NLO and LO results at 13 TeV for Ag defined as 

. ^ c^i\yt\ > \yi\) - (^{\yt\ < \yi\) on 

cT{\yt\>\yt\) + cy{\yt\<\yt\) 

are listed in table 3, which clearly demonstrates that NLO QCD effects cannot be neglected, 
once again, in the predictions of the asymmetries. For tiW^ and tiH production, an 
asymmetry is actually generated only at NLO. Furthermore, NLO QCD corrections change 
sign and increase by a factor ~ 7 the asymmetry in tiZ production and they decrease it 
by a factor larger than two in tty production. Thus, NLO results point to the necessity 
of reassessing the phenomenological impact of the tty signature, which is based on a LO 
calculation [59]. Moreover, we have also checked that for pr( 7 ) > 50GeV both the LO and 
NLO central values of the asymmetry are very similar (within 5 per cent) to the results in 
table 3, where pt( 7 ) > 20GeV. 

2.2 tiVV processes 

We start showing for all the tiVV processes the dependence of the NLO total cross sections, 
at 13 TeV, on the variation of the renormalisation and factorisation scales fir and yLj. This 
dependence is shown in fig. 11 and it is obtained by varying /r = = /rj by a factor eight 

around the central value yt = y,g (solid lines), y, = ya (dashed lines) and y = mt (dotted 
lines). Again, for all the processes and especially for those with a photon in the final state, 
we find that ya typically leads to smaller cross sections than yg and m*. For this class of 
processes we also investigated the effect of the independent variation of factorisation and 
renormalisation scales. We found that the condition y^ = yy captures the full dependence 
in the {yr,yf) plane in the range ^a/2 < yf,yr < 2^a- On the other hand, in the full 
ya/8 < yf,yr < 8ya region off-diagonal values might differ from the values spanned at 
yj — /ij.. 

In table 4 we list, for all the processes, LO and NLO cross sections together with PDF 
and scale uncertainties, and iL-factors for the central values. Again scale uncertainties 
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Figure 11. Cross sections of ttVV processes at 13 TeV. Comparison of NLO scale dependence in 
the interval ^'^/8 < fi < 8^'^ for the three different choices of the central value fj,g, Ha, rnt- 

are evaluated by varying independently the factorisation and the renormalisation scales in 
the interval fig/2 < < 2fig. The dependence of the LO and NLO cross sections on 

fi = fir = fJ-f shown in fig. 12 in the range fig/8 < fi < 8fig. As expected, for all the 
processes, the scale dependence is strongly reduced from LO to NLO predictions both in 
the standard interval fig/2 < fi < 2fig as well as in the full range fig/8 < fi < 8fig. For 
the central scale fi = fig, iC-factors are very close to unity. It is interesting to note that 
NLO curves display a plateau around fig/2 or fig/4:, corresponding to Ht/8 and Ht/^Q, 
respectively. 

We show now the impact of NLO QCD corrections for relevant distributions and we 
discuss their dependence on scale choice and its variation. For all the processes we have 
considered the distribution of the invariant mass of the top-quark pair and the pr and the 
rapidity of the (anti)top quark, of the top-quark pair and of the vector bosons. Again, given 
the large amount of distributions that is possible to consider for such a final state, we show 
only representative results. We remind the interested reader that additional distributions 
can be easily produced via the public code MadGraph5_aMC@NLO. 

For each figure, we display together the same type of distributions for the six different 
processes: M 77 , tiZZ, tiW~^W~, ttW^Z, tiW^^ and tiZ'^. We start with fig. 13, which 
shows the m{tt) distributions. The format of the plot is the same used for most of the 
distribution plots in subsection 2.1, where it is also described in detail. For m{tt) distri¬ 
butions, we notice features that are in general common to all the distributions and have 
already been addressed for ttV processes in subsection 2.1. For instance, the use oi fi = fia 
leads to NLO values compatible with, but systematically smaller than, those obtained with 
fi = fig. Conversely, the choice fi = mt leads to scale uncertainties bands that overlaps 


19 - 







13 TeV a[fb] 
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UW+W-[4f\ 



NLO 

o -1 1'v+3.8% +1.9% 
^■^^'-8.6% -1.8% 

1 1 Q/I“I“8-3% +2.3% 

±i.c5^ll 2% 2.4% 

10-261S 

+ 1.3% 
-1.3% 

LO 

0 io'7+36.1% +1.9% 

Z.±0^24 4% i.9% 

1 n '7Q“I“38.3% +2.2% 
±U.I0 25 4% 2.2% 

Q qqq+36.5% 
(5.5^55 24 5% 

+ 1.5% 
-1.6% 

AT-factor 

0.99 

1.10 

1.16 


13 TeV a[fb] 

tiw^z 

ttZj 


NLO 

A 1 r'^+9.8% +2.2% 

^.±O/_10 7% _i_Q% 

r '^'71+10.5% +1.8% 
0./ /I_i2.1% -1.9% 

6.734+(?;°^” 

+ 1.8% 
-1.4% 

LO 

o Q 01 +32,6% +2.3% 

0 .yzi _22 8% - 2 . 2 % 

c: nQn+38.0% +1.9% 

O.UOU-25 3% _i.9% 

6-i45+^22;f: 

+2.1% 

-2.0% 

AT-factor 

1.06 

1.14 

1.10 



Table 4. NLO and LO cross sections for ttVV processes at 13 TeV for p = The first uncertainty 
is given by the scale variation within /ig/2 < the second one by PDFs. The relative 

statistical integration error is equal or smaller than one permille. 



Figure 12. NLO and LO cross sections at 13 TeV. Scale dependence in the interval /r'^/8 < p, < 8/j,° 
with = /ig for the tiVV processes. 

with those obtained with fj, = fj,g. The NLO corrections in ttZZ production are very close 
to zero, for fi = fig, and very stable under scale variation (see also table 4). For all the 
processes, the two dynamical scales ^g and yield flatter iL-factors than those from the 
fixed scale mt- 

In fig. 14 we show the distributions for As for ttV processes (see fig. 4), these 

distributions receive large corrections in the tails. This effect is especially strong for the 
processes involving a photon in the final state, namely, ttjj, ttZ^ and Also, for all 


20 



















the three choices of /i employed here, ill-factors are not flat. Surprisingly, the i^-factors for 
tiZZ, ttW^Z and tiW~^W~ production show a larger dependence on the value of Pritt) 
when ^ is a dynamical quantity, as can be seen from a comparison of the first {p = Pg) and 
second (/i = pa) insets with the third insets {p = nit). From the fourth insets of all the six 
plots, it is possible to notice how the scale dependence at NLO for p = pg it is much larger 
than for p = pa- Exactly as we argued for ttV processes, NLO ttW+jets merged sample 
a la FxFx should be used for an accurate prediction of these tails. 

In fig. 15 we show the distributions for Prit)- Most of the features discussed for m{tt) 
in fig. 13 appear also for these distributions. The same applies to the distributions of the 
Pt of the two vector bosons, which are displayed in fig. 16. In the plots of fig. 16 and in 
all the remaining figures of this section we use the same format used in subsection 2.1 for 
fig. 10. Thus, differential iL-factors will not be explicitly shown. In the first and second 
inset we show the ratio of the distributions of the pr of the two vector bosons, respectively 
at NLO and LO accuracies. In the case of tt'j'j production, 71 is the hardest photon, while 
72 is the softest one. Similarly, in ttZZ production, Z\ is the hardest Z boson, while Z 2 is 
the softest one. As can be noticed, for each process this ratio is the same at LO and NLO 
accuracy and thus it is not sensitive to NLO QCD corrections. 

In fig. 17 we show the distributions for y{t) and y(t). The ttVV processes, with the 
exception of ttW~^W~ at LO exhibit a central asymmetry for top (anti-)quarks. Top 
quarks are more centrally distributed than top antiquarks in ttjj, ttW^j and tiZ'j produc¬ 
tions, while they are more peripherally distributed in tiZZ and ttW^Z production. In all 
the tiVV processes, NLO QCD corrections lead to a relatively more peripheral distribution 
of top quarks than antiquarks. This effects yield to a non-vanishing central asymmetry 
for tiW~^W~ production and almost cancel the LO central asymmetry of tiZ^ production. 
Here, we refrain to present results for the central asymmetries of tiVV processes, since it 
is extremely unlikely that at the LHC it will be possible to accumulate enough statistics to 
perform these measurements. 

In fig. 18 we show the distributions for y(Li) and y(V 2 ). Comparing the first and second 
insets, only small differences can be seen for the ratios of the distributions at LO and NLO. 
Thus, unlike for the top quark and antiquark, the rapidity of the first and the second vector 
boson receive NLO relative differential corrections that are very similar in size. Both in the 
distributions of the rapidities of the top (anti)quark and of the vector bosons, NLO QCD 
corrections in general induce non-flat iL-factors, also with the use of dynamical scales.® 


^Analytically, this process is supposed to give an asymmetry. Numerically, it turns out that it can be 
safely considered as zero. 

®We explicitly verified it and it can be easily reproduced via the public version of Mad- 
GrAPH5_aMC@NLO, which has also been used for the phenomenological study presented here. 
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Figure 13. Differential distributions for the invariant mass of top-quark pair, The format 

of the plots is described in detail in subsection 2.1. 
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Figure 14. Differential distributions for the Pt of top-quark pair, pT{tt). The format of the plots 
is described in detail in subsection 2.1. 
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Figure 15. Differential distributions for the of top-quark, Prit). The format of the plots 
described in detail in subsection 2.1. 












































































































































Figure 16. Differential distributions for the of the first and second vector boson, PriVi) and 
Pt{V2). 
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Figure 17. Differential distributions for the rapidity of the top quark and antiquark, y{t) and y{t). 
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Figure 18. Differential distributions for the rapidity of the first and second vector boson, y{Vi) 
and y{V2). 
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Figure 19. NLO and LO cross sections for tttt production at 13 TeV. Comparison of the scale 
dependence in the interval ^'^/8 < /r < for the four different choices of the central value /i°: fig, 

fia, fiY^, 2mt. 

2.3 tttt production 

In this section we present results for tttt production. We start by showing in fig. 19 the 
scale dependence of the LO (blue lines) and NLO (red lines) total cross section at 13 TeV. 
As for the previous cases, we vary fi = fi^ = fif hy a factor eight around the central value 
fi = fig (solid lines), fi = fia (dashes lines) and, due to the much heavier final state, fi = 2mt 
(dotted lines). In this case we also show with a dot-dashed line the dependence of the NLO 
cross section on an alternative definition of average scale fi]^^ = where 

possible additional partons appearing in the final state do not contribute. 

As expected, predictions relative to fig and fi^^ are very close. Conversely, fia and 

show a non-negligible difference. Note that the value of fia and fi}^^ is the same for 
Born and and virtual contributions for any kinematic configuration. Thus, the difference 
between dashed and dot-dashed lines is formally an NNLO effect that arise from differences 
in the scale renormalisation for real radiation events only. To investigate the origin of this 
effect, we have explicitly checked that the difference is mainly induced by the corresponding 
change in the renormalisation scale and not of the factorisation scale. Similar behaviour is 
also found in tiV and tiVV processes, yet since the masses of the final-state particles are 
different and the as coupling order lower, fig and lines are more distant than in tttt 
production. 

Since the LO cross section is of 0(ag), it strongly depends on the value of the renor¬ 
malisation scale, as can be seen in fig. 19. This dependence is considerably reduced at NLO 
QCD accuracy in the standard interval fig/2 < fi < 2fig. Conversely, for fi < fig/4: the value 
of the cross section falls down rapidly, reaching zero for fi ~ lig/8. This is a signal that 
in this region the dependence of the cross section on /x is not under control. Qualitatively 
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similar considerations apply also for the different choices of scales, as can be seen in fig. 19. 

In eqs. (2.4) and (2.5), we list the NLO and LO cross sections evaluated at the scale 
IJ, = together with scale and PDF uncertainties. As done in previous subsections, scale 
uncertainties are evaluated by varying the factorisation and renormalisation scales in the 
standard interval fJ-g/2 < As a result the total cross section at LHC 13 TeV 

for the fj, = fj,g central scale choice reads 

^NLO = 13.31+25.8%+5.8% fb, (2.4) 

ctlo = WM+imttS fb, (2.5) 

iP—factor = 1.22. (2-6) 

Different choices for the central value and functional form of the scales, as well as the 
interval of variation, lead to predictions that are compatible with the result above, see also 
e.g. [16]. 

We now discuss the effect of NLO QCD corrections on differential distributions. We 
analysed the distribution of the invariant mass, the pr and the rapidity of top (anti-)quark 
and the possible top-quark pairs. Again, given the large amount of distributions, we show 
only representative results. All the distributions considered and additional ones can be pro¬ 
duced via the public code MadGraph5_aMC@NLO. For this process the scale dependence 
of many distributions has been studied also in [16] and our results are in agreement with 
those therein. In fig. 20 we show plots with the same formats as those used and described 
in the previous sections. Specifically, we display the distributions for the total pj. of the two 
hardest top quark and antiquark (Pt(HH))) their invariant mass the rapidity of 

the hardest top quark y{ti) and the invariant mass of the tttt system Also, in 

the last plot of fig. 20, we show the Pt distributions of the hardest together with the softest 
top quarks, priti) and Priti)) and their ratios at NLO and LO. 

We avoid repeating once again the general features that have already been pointed out 
several times in the previous two sections; they are still valid for ttti production. Here, we 
have found, interestingly, that NLO corrections give a sizeable enhancement in the threshold 
region for It is worth to notice that also for this process NLO QCD corrections 

are very large in the tail of the Pt(HH) distribution, especially with the use of dynamical 
scales. We have verified that in these regions of phase space the qg —)• ttttq contributions 
are important. Finally, as can be seen in the last plot, we find that the ratios of Pt(H) and 
Py(t 2 ) distributions are not sensitive to NLO QCD corrections. 
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Figure 20. Differential distributions for tttt production. 
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2.4 Total cross sections from 8 to 100 TeV 

In addition to the studies performed for the LHC at 13 TeV, in this subsection we discuss 
and show results for the dependence of the total cross section on the energy of the proton- 
proton collision. In fig. 21 NLO QCD total cross sections are plotted from 8 to 100 TeV, 
as bands including scale and PDF uncertainties. The corresponding numerical values are 
listed in table 5. As usual, central values refers to /i = fig, and scale uncertainties are 
obtained by varying independently and /ij in the standard interval ^g/2 < fif, < 2//^. 

In the upper plot of fig. 21 we show the results for tiH production and tiV processes, 
whereas tttt production and ttVV processes results are displayed in the lower plot. In both 
plots we show the dependence of the iL-factors at = ^g on the energy (the first and the 
second inset). The first insets refer to processes with zero-total-charge final states, whereas 
the second insets refer to processes with charged final states. The very different qualitative 
behaviours between the two classes of processes is due to the fact that the former include 
already at LO an initial state with gluons, whereas the latter do not. The gluon appears 
in the partonic initial states of charged processes only at NLO via the (anti)quark-gluon 
channel. At small Bjorken-x’s, the gluon PDF grows much faster than the (anti)quark 
PDF. Thus, increasing the energy of the collider, the relative corrections induced by the 
(anti)quark-gluon initial states leads to the growth of the iP-factors and dominates in their 
energy dependence. Also, as can be seen in fig. 21 and table 5, these processes present a 
larger dependence on the scale variation than the uncharged processes. 

The differences in the slopes of the curves in the main panels of the plots are also 
mostly due to the gluon PDF. Charged processes do not originate from the gluon-gluon 
initial state neither at LO nor at NLO. For this reason, their growth with the increasing of 
the energy is smaller than for the uncharged processes. All these arguments point to the 
fact that, at 100 TeV collider, it will be crucial to have NNLO QCD corrections for ttW^, 
ttW^j and ttW^Z processes, if precise measurements to be compared with theory will be 
available. 

The fact that tttt production is the process with the rapidest growth is again due 
to percentage content of gluon-gluon-initiated channels, which is higher than for all the 
other processes, see fig. 22. From the left plot of fig. 21, it is easy also to note that the 
scale uncertainty of tttt production is larger than for the tiVV processes. In this case, 
the difference originates from the different powers of as at LO; tttt production is of O(a^) 
whereas ttVV processes are of ©(a^a^). 
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Figure 21. NLO total cross sections from 8 to 100 TeV. The error bands include scale and PDF 
uncertainties (added linearly). The upper plot refers to tiV processes and tiH production, the lower 
plot to ttVV processes and titt production. For final states with photons the Pt(,i) > 20 GeV cut 
is applied. 
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Figure 22. Relative contribution of the gg channel to the total cross section at LO for 
tiV, tiH, tiVV and tttt processes for pp collisions from 8 to 100 TeV centre-of-mass energy. For 
final states with photons the Pril) > 20 GeV cut is applied. 
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a [pb] 
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13 TeV 

14 TeV 

25 TeV 
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U.ZZD-11 9% _3 
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n 'tqo + 12.7% +2.1% 
^■'°^-13.5% -2.4% 

9 +1.6% 
-1.9% 

0 9/5+14.2% +1.6% 
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11.1 1-12 7% -1.4% 

9n 84+14.9% +1.1% 
^'J.0+_^2.5% -1.3% 
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Table 5. NLO cross sections for tiVV,tttt,ttV,tiH processes using the geometrical average scale. The first uncertainty is given by scale variation, 
the second by PDFs. For final states with photons the Pt{i) > 20 GeV cut is applied. The cross sections for the four final-state particle processes 
are calculated with percent accuracy, whereas for the processes with three final-state particles with per mill. 



3 Analyses of ttH signatures 

In this section we provide numerical results for the contributions of signal and irreducible 
background processes to two different classes of ttH signatures at the LHC. In subsection 
3.1 we consider a signature involving two isolated photons emerging from the decay of the 
Higgs boson into photons, H —?■ 77 . In subsection 3.2 we analyse three different signatures 
involving two or more leptons, where tiH production can contribute via the H —?■ ZZ*, 
H —)• WW* and H —)■ decays. We perform both the analyses at 13 TeV and we adopt 

the cuts of [27].® The preselection cuts, which are common for both the analyses, are: 

PT(e)>7GeV, |? 7 (e) < 2.5| , pT(/u)>5GeV, |r/(^)|<2.4, 

17 ( 7)1 <2.5, p^(j)>25GeV, | 7 (j)|< 2 . 4 , (3.1) 

where jets are clustered via anti-fer algorithm [60] with the distance parameter R = 0.5. 
Event by event, only particles satisfying the preselection cuts in eq.(3.1) are considered 
and, for each jet j and lepton if AR(j,i) < 0.5 the lepton i is clustered into the jet j. 
With the symbol £, unless otherwise specihed, we always refer to electrons(positrons) and 
(anti)muons, not to r (anti)leptons. 

All the simulations for the signal and the background processes have been performed 
at NLO QCD accuracy matched with parton shower effects (NLO + PS). Events are gener¬ 
ated via MadGraph5_aMC@NLO, parton shower and hadronization effects are realised in 
PythiaS [39], and jets are clustered via Fast Jet [61].^® Unless differently specihed, decays 
of the heavy states, including r leptons, are performed in PythiaS. In the showering, only 
QGD effects have been included; QED and purely weak effects are not included. Further¬ 
more, multi-parton interaction and underlying event effects are not taken into account. 

In order to discuss NLO effects at the analysis level, in the following we will also report 
results for events generated at LO accuracy including shower and hadronization effects 
(LO -I- PS). As done for the hxed-order studies in section 2, LO J- PS and NLO + PS 
central values are evaluated at fj,f = = fJ-g and scale uncertainties are obtained by 

varying independently the factorisation and the renormalisation scale in the interval fj,gj2 < 

< 2//g. 

3.1 Signature with two photons 

The present analysis focuses on the Higgs boson decaying into two photons in ttH produc¬ 
tion, which presents as irreducible background the ttjj production. In our simulation, top 
quark pairs are decayed via Madspin for both the signal and the background, whereas the 
loop-induced 77 —)• 77 decay is forced in PythiaS and event weights are rescaled by the 
branching ratio BR(i7 —)> 77 ) = 2.28 x 10“^, which is taken from [62]. 

In this analysis, at least two jets are required and one of them has to be 6 -tagged. In 
addition, the following cuts are applied: 

100 GeV < m( 7 i 72 ) < 180 GeV, Pt{ji) > ^ Pt{ 12 ) > 25GeV, 

®In our simulation we do not take into account particle identification efficiencies and possible misidenti- 
fication effects. 

^°In our simulation, fo-tagging is performed by looking directly at B hadrons, which we keep stable. 
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13 TeV cr[fb] 

ttH X BR(iL —¥ 77 ) tt'y'y 


ttH(H 77 ) £t 77 

NLO 

1 191+6'09'o + 2 . 1 % 1 466+®-'^^" 

-^•-^^-^-9,4% - 2 . 6 % - 1 . 8 % 

NLO+PS 

0.194+=;®^° ± 0.002 0-374+“'jP^ ± 0.004 

LO 

, ^5,7+35.5% +2.0% 1 o-in+37.0% +1.7% 

-24.2% -2.1% ^•'^^'^-24.8% -1.8% 

LO+PS 

0.i72t^®;(P° +2:2% ± 0.001 0.310+1®;*! ± 

K 

1.10 1.09 


1.13 ±0.01 1.21 + 0.01 


Table 6. NLO and LO cross sections for ttH{H —> processes at 13 TeV. The first 

uncertainty is given by scale variation, the second by PDFs. The assigned error is the statistical 
Monte Carlo uncertainty. 


Ai?( 7 i, 72 ), AR{'yi^ 2 ,j) > 0.4, AR{'yi^ 2 j) > 0.4 , prih) > 20GeV, (3.2) 


and an additional cut 


AR{iiJj)>0A (3.3) 

is applied if leptons are more than one. With 71 and 72 we respectively denote the hard 
and the soft photon, analogously £i indicates the hardest lepton. Cuts on lepton(s) imply 
that the fully and semileptonic decays of the top-quark pair are selected. 

Results at LO -I- PS and NLO -|- PS accuracy are listed in table 6 for the signal and the 
tt'j'y background. Also, we display hxed order results (LO, NLO) at production level only, 
without including top decays, shower and hadronization effects. In order to be as close as 
possible to the analyses level, we apply the cuts in eqs. (3.1) and (3.2) that involve only 
photons. Thus, the difference between LO and NLO results of ttj'j in tables 4 and 6 are 
solely due to these cuts. 

In table 6 , we show global iL-factors both at hxed order {K := NLO/LO) and in¬ 
cluding decays, shower and hadronization effects, and all the cuts employed in the analysis 
:= NLO-|-PS/LO-I-PS). Comparing and K it is possible to directly quan¬ 
tify the difference between a complete NLO simulation (K^^) and the simulation typically 
performed at experimental level, i.e., a LO-|-PS simulation rescaled by a iP-factor from 
production only {K). As shown in table 6 , e.g., the second approach would underestimate 
the prediction for ttyy production wrt. a complete NLO -|- PS simulation. This difference is 
not of particular relevance at the level of discovery, which mostly relies on an identihcation 
of a peak in the m( 7 i 72 ) (see also hg. 23), but could be important in the determination of 
signal rates and in the extraction of Higgs couplings. Conversely, the difference between K 
and is much smaller for the signal. 

In hg. 23 we show representative differential distributions at NLO -|- PS accuracy for the 
signal (red) and background (black) processes. In the two insets we display the differential 
iL-factor for the signal and the background (iL^®^) using the same layout and 

conventions adopted in the plots of section 2. In particular, we plot the invariant mass of 
the two photons ( 771 ( 7172 )) their distance (Ai 2 ( 7 i, 72 )) and the transverse momentum of the 
hard (Pt(7i)) and the soft (Pt( 72 )) photon. We note that predictions for key discriminating 
observables, such as the Ai?( 7 i, 72 ) and Pt( 72 ) are in good theoretical control. 
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Figure 23. Differential distributions for signal and background processes for the diphoton analysis. 

3.2 Signatures with leptons 

This analysis involves three different signatnres and signal regions that inclndes two or 
more leptons and it is specifically designed for ttH prodnction with snbseqnent H —)• ZZ*, 
H —)• WW* and H —)> decays. In the simnlation, all the decays of the massive 

particles are performed in PythiaS. In the case of the signal processes, the Higgs boson 
is forced to decay to the specific final state {H —)• ZZ*, H —)• WW* or H ^ and 

event weights are rescaled by the corresponding branching ratios, which are taken from 
[62]: BR{H WW*) = 2.15 x 10-\ BR{H ZZ*) = 2.64 x IQ-^, BR{H r+r”) = 
6.32 X 10“^. The isolation of leptons from the hadronic activity is performed by directly 
selecting only prompt leptons in the analyses, i.e., only leptons emerging from Z, W or 
from r leptons which emerge from Z, W or Higgs bosons. 

We consider as irredncible backgronnd the contribntion from ttW^, tiZ/j*, ttW~^W~, 

observed that applying hadronic isoiation cuts as done in [27] we obtain resuits with at most 
10% difference with those presented here by seiecting prompt ieptons. il-factors are independent of the 
appiication of hadronic isoiation cuts. 
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tiZZ^ ttW^Z and tttt production. Precisely, with the notation ttZ/'j* we mean the full 
process = e,/r,r), where Z and photon propagators, from which the pair 

emerges, can both go off-shell and interfere.All the processes, with the exception of 
ttZ/'j*, have been also studied at fixed-order accuracy in section 2. 

In the analyses the following common cuts are applied in order to select at least two 
leptons 

mihh) > 12, AR{ii,£j) > 0.4. (3.4) 

Then, the three signatures and the corresponding signal regions are defined as described in 
the following: 

• Signal region one (SRI): two same-sign leptons 

Exactly two same-sign leptons with with pt{^) > 20 GeV are requested. The event is 
selected if it includes at least four jets with one or more of them that are 6-tagged. Further¬ 
more it is required that pt(^i) +Pt(^ 2 ) + > 100GeV and, for the dielectron events, 

|m(e^e^) — mz\ > 10 GeV and El™®® > 30 GeV, in order to suppress background from 
electron sign misidentification in Z boson decays. 

• Signal region two (SR2): three leptons 

Exactly three leptons with pt i^i) > 20 GeV, pt{(- 2 ) > 10 GeV, pt{(-^ = e(//)) > 7(5) GeV 
are requested. The event is selected if it includes at least two jets with one or more of them 
that are 6-tagged. For a Z boson background suppression, events with an opposite-sign 
same-flavour lepton pair are required to have — mz\ > 10 GeV. Also, for this kind 

of events if the number of jets is equal or less than three, the cut El™'®® > 80 GeV is applied. 

• Signal region three (SR3): four leptons 

Exactly four leptons with pt{£i) > 20 GeV, pt(^ 2 ) > 10 GeV, pr(^ 3,4 = e(/i)) > 7(5) GeV 
are requested. The event is selected if it includes at least two jets with one or more of 
them that are 6-tagged. Also here, for a Z boson background suppression, events with an 
opposite-sign same-flavor lepton pair are required to have |m(7+£“) — mz\ > 10 GeV. 

For both signal and background processes, results at LO -|- PS and NLO -|- PS accuracy 
as well as Ef^^-factors are listed in table 7 for the three signal regions. Also, for each process 
we display the value of the global AT-factor (listed also in section 2), which does not take 
into account shower effects, cuts and decays. A posteriori, we observe that in these analyses 
the AT-factors are almost insensitive of shower effects and the applied cuts. This is evident 
from a comparison of the values of K and in table 7, where the largest discrepancy 
stems from the tiZ /^* process in SRI. We also verified, with the help of Madspin, that 
results in the SR3 (SR2 for tiW ^) do not change when spin-correlation effects are taken 
into account in the decays.It is important to note that, a priori, with different cuts 

principle also ttW^ and tiZ') production can contribute to the signatures specified in the following. 
However, they are a small fraction of ttW and ttZ production and indeed are not taken into account in the 
analyses of [27]. 

^®To this purpose, we excluded Higgs boson propagators in order to avoid a double count of the ttH{H —> 
contributions. 

^^SR2 and especially SRI involves a rich combinatoric of leptonic and hadronic Z, W and r decays, which 
render the simulation with spin-correlation non-trivial. However, we checked also here for representative 
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0.0087i1^!|% 1+!: ± 0.00007 

K = 1.22 

^PS 

1.26 ± 0.02 

1.27 ± 0.02 

1.26 ± 0.03 

13 TeV cr[ah] 


SRI 

SR2 

SR3 


NLO+PS 

9 60±^'^^ +1.8% m A Ag 

5 02 + ^'^^ +1-8% + 0 04 

0.2491!“% 1!“% ± 0.009 

tiZZ 

LO+PS 

g „4 + 36.3% +1.9% , Q Q2 
-24.5% -1.9%^^'^ 

c r)Q + 35.9% +1.9% 1 

5-08-24 3% -1.9% 

0.25o1“!|% +“!: ± 0.004 

K - 0.99 

^PS 

0.99 ± 0.01 

0.99 ± 0.01 

1.00 ± 0.04 


NLO+PS 

g2 a+ 9,0% +2.2% , g 7 

‘^^•^-10.2% -1.6 %*^'^ 

27 g + 9.2% +2.3% ,05 

^+^-10.3% -1.7%^^’^ 

0-91I!/!: 1+!: ± 0.02 

ttw^z 

LO+PS 

gg 2 + 32-2% +2.4% , g 3 
^^•^-22.6% -2.3% ^ 

2g 4+32.0% +2.4% ,02 
^5-4-22 5% -2.2% 

0.8931“!“% tHl ± 0.009 

K - 1.06 

^PS 

1.03 ± 0.01 

1.06 ± 0.02 

1.02 ± 0.02 


Table 7. NLO and LO cross sections for signal and background processes for tiH to multileptons 
at 13 TeV. The first uncertainty is given by scale variation, the second by PDFs. The assigned error 
is the statistical Monte Carlo uncertainty. 

and/or at different energies, K and could be in principle different and spin correlation 
effects may be not negligible. Thus, a genuine NLO+PS simulation is always preferable. 


cases that spin-correlation effects do not sensitively alter the results. 
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4 Conclusions 


In this paper we have presented a thorough study at NLO QCD accuracy for tiV and tiVV 
processes as well as for ttH and tttt production within the same computational framework 
and using the same input parameters. In the case of ttVV processes, with the exception of 
ft 77 production, NLO cross sections have been studied for the hrst time here. Moreover, 
we have performed a complete analysis with realistic selection cuts on hnal states at NLO 
QCD accuracy including the matching to parton shower and decays, for both signal and 
background processes relevant for searches at the LHC for the ttH production. Specihcally, 
we have considered the cases where the Higgs boson decays either into leptons, where tiV 
and ttVV processes and titi production provide backgrounds, or into two photons giving 
the same signature as ti'j'j production. 

We have investigated the behaviour of hxed order NLO QCD corrections for several dis¬ 
tributions and we have analysed their dependence on (the dehnition of) the renormalisation 
and factorisation scales. We have found that QCD corrections on key distributions cannot 
be described by overall iL-factors. However, dynamical scales in general, even though not 
always, reduce the dependence of the corrections on kinematic variables and thus lead to 
flatter iL-factors. In addition, our study shows that while it is not possible to identify 
a “best scale” choice for all processes and/or differential distributions in tiV and ttVV, 
such processes present similar features and can be studied together. For all the processes 
considered, NLO QCD corrections are in general necessary in order to provide precise and 
reliable predictions at the LHC. In particular cases they are also essential for a realistic 
phenomenological description. Notable examples discussed in the text are, e.g., the giant 
corrections in the tails of Priti^ distributions for ttV processes and the large decrement of 
the top-quark central asymmetry for ti'j production. In the case of future (hadron) colliders 
also inclusive cross sections receive sizeable corrections, which lead, e.g., to iL-factors larger 
than two at 100 TeV for ttV and tiVV processes with a charged hnal state. 

In the searches at the LHC for the ttH production with the Higgs boson decaying either 
into leptons or photons, NLO QCD corrections are important for precise predictions of the 
signal and the background. We have explicitly studied the sensitivity of NLO+PS QCD 
corrections on experimental cuts by comparing genuine NLO+PS QCD predictions with 
LO+PS predictions rescaled by global iP-factors from the hxed order calculations without 
cuts. A posteriori, we have verihed that these two approximations give compatible results 
for analyses at the 13 TeV Run-H of the LHC with the cuts specihed in the text. A priori, 
this feature is not guaranteed for analyses with different cuts and/or at different energies. 
In general, a complete NLO+PS prediction for both signal and background processes is 
more reliable an thus preferable for any kind of simulation. 

All the results presented in this paper have been obtained automatically in the publicly 
available MadGraph5_aMC@NLO framework and they can be reproduced starting from 
the input parameters specihed in the text. 
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